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Abstract: We revisit the prospects for detecting the Sunyaev Zel'dovich (SZ) effect in- 
duced by dark matter (DM) annihilation or decay. We show that with standard (or even 
extreme) assumptions for DM properties, the optical depth associated with relativistic 
electrons injected from DM annihilation or decay is much smaller than that associated 
with thermal electrons, when averaged over the angular resolution of current and future 
experiments. For example, we find: tdm ~ 10 -9 — 10 -5 (depending on the assumptions) 
for m x = 1 GeV and a density profile p oc r _1 for a template cluster located at 50 Mpc 
and observed within an angular resolution of 10", compared to r t h ~ 10 -3 — 10~ 2 . This, 
together with a full spectral analysis, enables us to demonstrate that, for a template clus- 
ter with generic properties, the SZ effect due to DM annihilation or decay is far below 
the sensitivity of the Planck satellite. This is at variance with previous claims regarding 
heavier annihilating DM particles. Should DM be made of lighter particles, the current 
constraints from 511 keV observations on the annihilation cross section or decay rate still 
prevent a detectable SZ effect. Finally, we show that spatial diffusion sets a core of a few 
kpc in the electron distribution, even for very cuspy DM profiles, such that improving the 
angular resolution of the instrument, e.g. with ALMA, does not necessarily improve the 
detection potential. We provide useful analytical formulae parameterized in terms of the 
DM mass, decay rate or annihilation cross section and DM halo features, that allow quick 
estimates of the SZ effect induced by any given candidate and any DM halo profile. 
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1. Introduction 

The enigma of the origin of dark matter (DM) is a longstanding issue. Many experi- 
mental strategies in particle physics, astrophysics and cosmology have been developed to 
confirm/infirm the existing scenarios. One appealing possibility, strongly motivated by 
independent issues in particle physics beyond the standard model (BSM), is that DM is 
made of weakly interacting particles (WIMPs). The most popular models, like the super- 
symmetric (SUSY) neutralino (see e.g. |IL |2| for reviews), have the interesting property 
of self-annihilation to ordinary matter, which provides windows for the detection of the 
annihilation products. There is actually an impressive wealth of WIMP models, from the 
GeV-TeV mass-scale as in the SUSY paradigm down to the sub-GeV mass-scale, the latter 
class allowing, in particular, to solve the so-called cusp and subhalo problems (e.g. §). 
Since the DM annihilation or decay rate increases with the DM number density, the best 
targets are expected to be the centers of (sub)galaxies or galaxy clusters. This quest for 
annihilation or decay signals makes use of different messengers (e.g. 7-rays, antimatter 
cosmic rays and neutrinos): it is referred to as indirect detection of DM (for a pedagogical 
review, see Q). There are of course other ways to learn about the microscopic properties of 
DM, for instance by detecting direct interactions of DM in ground based experiments Q, 
or by producing DM directly in particle colliders. All these direct and indirect detection 
methods have been studied at length in the literature. 
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Yet, a few years ago, it was pointed out in a series of papers ||, |(| 0] a new comple- 
mentary method to look for DM, namely the Sunyaev-Zel'dovich (SZ) effect ]|, [|. Indeed, 
DM annihilation or decay in galaxy clusters may inject relativistic electrons and positrons 1 
that should experience inverse Compton scattering with the CMB photons and, therefore, 
if numerous enough, generate a deviation to the blackbody spectrum. The question that 
we address in this paper is whether or not this deviation is observable in light of the con- 
tradictory results found by pL pi m on the one hand, and by [10] on the other hand. In 
particular, we investigate whether the Planck satellite [11], which was recently launched, 
and the ALMA facilities Jl^, [l3| will have the sensitivity to detect the SZ signal induced 
by DM annihilation or decay and could constrain the DM properties. 

Our approach aims at providing more analytical insights than [1C] on the crucial im- 
pact of the experimental angular resolution on the DM-induced SZ predictions, on the 
spectral analysis itself and on spatial diffusion effects. We also consider the SZ effect due 
to the thermal electrons lying in clusters and observed in X-ray measurements, which has 
been studied for a long time (see [14] for a review) and which constitutes an important 
background to such an exotic signature. To proceed, we will consider a template galaxy 
cluster located at 50 Mpc, with properties representative of all nearby clusters. Our dis- 
cussion, that we will keep as generic as possible, is organized as follows. In Sec. we will 
recall the principle of the SZ calculation. In Sec. || we will derive the optical depths for 
both DM-induced and thermal electrons in an exact numerical manner and, likewise, pro- 
vide analytical formulae to allow quick calculation of the SZ effect for any cluster and DM 
properties. In Sec. we will then use these results together with a full spectral analysis 
to compute the signal due to the DM-induced electrons at the transition frequency where 
the thermal component is negligible, and show that the DM contribution to the SZ will 
be quite hardly observable with current and even future radio experiments. In Sec. ||, we 
will close the case by considering spatial diffusion effects, showing analytically that they 
induce a core in the electron distribution whatever the cuspiness of the DM profile. We 
will conclude in Sec. 0. 



2. Principle 



There are basically two main approaches to calculate the thermal or relativistic SZ effect, 



the radiative transfer method proposed in e.g. [15] and early formalized in [16], and the 
covariant Boltzmann equation formalism [17, 18, 19, 20, 21]. Recently, however, these 



approaches have been demonstrated to be equivalent [22, 23 1, which settles a self-consistent 
and unique framework for the SZ calculations. In the single scattering approximation, the 
distortion of the CMB intensity / at a photon energy E k is given by the following equation: 



AIJE k ) = t 



dt (P^t) - 5(t - 1)) I°(E k /t) 



(2.1) 



where I® is the undistorted blackbody intensity, Pi is the frequency redistribution function 
for a single electron-photon interaction and t is the ratio of scattered to unscattered photon 



1 Hereafter, we will use electron for electron or positron, indifferently. 
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energy. The integrand encodes the spectral features: the first term describes the spectral 
part that is shifted from energy E^/t to E^, while the second term describes the part that 
is removed from energy Ek to others, r, which weights the amplitude of the distortion, 
is the optical depth which is related to the electron density in the cluster through the 
line-of-sight integral 

r = a T J dln e (x) . (2.2) 

Here ctt is the non-relativistic Thomson cross section that fully characterizes interactions 
of CMB photons with electrons of Lorentz factor 7 e < 10 8 . As far as r < 1, the single 
scattering approximation is fully justified [23|. 



In order to estimate the potential of a radio experiments to detect SZ signals, however, 
one needs to average the optical depth within the angular resolution of the apparatus, 



<7T 

AO(/i res ) 7AQ(/i ros ) 



< r Vr CS = A o/. ^ / dn I dln e (x) , (2.3) 



where /i res = cos(# res ) characterizes the angular resolution 6 Tes . 

Disregarding for the moment the spectral aspects, the observation of a relativistic SZ 
signal on top of the thermal contribution would roughly imply that the optical depth of the 
additional electron component is sizable compared to the one associated with the thermal 



electrons. A mere comparison of the electron densities in the cluster center, as done in [10], 
is indicative at th order. Nevertheless, because of the angular averaging, this could be 
misleading. Indeed, the SZ signal is in any case smeared down by angular resolution effects. 

In Sec. ||, we will study the effects of this averaging procedure, which is actually very 
common in studies on the indirect detection of DM with 7-rays (e.g. |24|] ). We will consider 
angular resolutions between ~ 0.1" and 1', noticing that Planck pi and ALMA fll2| , [l3|| 
could reach resolutions of ~ 5' and ~ 1", respectively. 

Throughout the paper, moreover, we will consider the spatial distribution of the ther- 
mal electrons to be a spherical cored isothermal: 



„o 

n e,th 



n e ,th(r) = ^ . ( 2 - 4 ) 

1 + 



where we will fix nP e th to 0.01 cm -3 , and r c to 400 kpc. These are typical values for nearby 
clusters that will partly characterize our template cluster. 

For completeness, it is worth recalling that the SZ effect is featured by the shift of 
the low frequency part of the CMB spectrum towards higher frequencies, but much more 
marginally by that of the high frequency part towards lower frequencies. This asymmetry 
drives AI to be negative at low frequencies and positive at high frequencies. The transition 
frequency, at which AI = 0, is intimately related to the energy spectrum of the involved 
electron gas. Because the thermal and DM-induced components are found in different 
energy ranges, the associated transition frequencies are different. Therefore, not only do 
we need to compare the optical depths of both contributions, that weight the amplitudes 
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of the corresponding distortions, but also to study the spectral features of AI. For the 
spectral study, we will mainly focus on AI X at the transition frequency of the thermal 
component (A/ t h = 0) in Sec. |j. 

3. Injected electron density from DM annihilation or decay and optical 
depths 

The injection rate of relativistic electrons (and positrons) from DM annihilation or decay 
characterizes a nominal source term 

Q na {E,x) = N a n (Pl^lY T(E) , (3.1) 

where iVo is the total number of electrons and positrons injected per annihilation or decay 
in the relevant energy range, J-(E) is the energy spectrum normalized to unity, and m x is 
the DM particle mass. The index n is 1 or 2 for DM decay or annihilation, respectively. 
Therefore, a\ = T x is a decay rate, while ct2 = 8{av)/2 (5 = 1 or 1/2 for Majorana or 
Dirac fermions, 1 for bosons) is an annihilation rate. The DM mass density profile p 7 is 
indexed by its inner logarithmic slope 7, and is usually written like a spherical symmetric 
component with a scale radius r s and a scale density p s [^] 

pJr) = Ps(rs/r) ^ = Ps Ur) . (3.2) 

The so-called Navarro-Frenk- White (NFW) profile [^] corresponds to (a, (3, 7) = (1,3,1), 
while the Moore profile [^7| is even more cuspy with (1.5,3,1.5). Notice that a cored 
isothermal profile with (2, 2, 0) is equivalent to the thermal electron distribution of Eq. ($2.4), 
with r s = r c ; this may be kept in mind with further benefit. We will study all of the three 
aforementioned profiles in the following, for both annihilating and decaying DM models. 
The typical values of the scale parameters found for galaxy clusters with masses ~ 1O 15 M 
in cosmological N-body simulations are r s ~ 400 kpc and p s ~ 0.05GeV.cm~ 3 (e.g. pq] ). 
We will adopt these parameters for the template cluster that we will use for our calculations 
throughout this paper. We summarize them in Tab. ||. 

distance i? v ; r r s and r c p s thermal e~ density 
(Mpc) (kpc) (kpc) (GeV/cm 3 ) (cm" 3 ) 
50 2000 400 O05 O01 

Table 1: Properties of the template cluster used in this paper. 

After their injection in the intracluster medium, relativistic electrons and positrons 
diffuse in space and momentum. The main processes that come into play are the energy 
losses, with a typical timescale ~ 300 Myrs, and the spatial diffusion due to the scattering 
on the magnetic inhomogeneities. It was shown in previous analyses that since the relevant 
spatial diffusion scale is ~ kpc, which is just a bit larger than typical resolution scales, 
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one can neglect spatial diffusion at first order |ffl]. We will stick to this approximation, for 
which, in steady state, the diffusion equation reduces to 



dE \ v ' dE 

where b{E) = —dE/dt is the energy loss rate. This equation is easily solved: 



N Q a n ( p^(x) 



, , , dE s T(E s 

b(E) \m x J J E V 

Now, to compute the optical depth, as defined in Eq. ( |2.2[ ), we need to integrate this 
differential electron density over energy. Using the previous equation, this gives: 



n n m^- 



LI mL m) - <3 - 5) 

We can further express the energy loss rate as b(E) = (Eq / 't\ oss ) / 'g(E) , where Eq = 1 
GeV, 

"Hoss is the typical energy loss timescale, and g(E) is a dimensionless function that 
encodes the energy dependence of the energy loss rate. The previous equation is then more 
conveniently rewritten as 

n n ^(x) = N a n n rioss ? , (3.6) 

where we have defined 

T = ^g(E) [~ dE s T{E s ) . (3.7) 

J Emm ^0 JE 

Note that T < m x /Eo for Coulomb losses (g(E) oc est), while T < Eo/E m i n for inverse 
Compton losses (g(E) oc (E/ Eq)~ 2 ), which is in any case < 10 3 . 

Now, to go further in our calculation of the optical depth, we need to compute the 
average density of electrons within the angular resolution of the telescope. 

3.1 Line-of-sight averaged optical depth and electron density 

As stated in Sec. || through Eq. the optical depth must be averaged within the 

experimental angular resolution. The full line-of-sight integral is not analytical in most of 
cases. We can reexpress Eq. ( |2.3| ) in such a way that an averaged electron density explicitly 
appears: 

(r^)=2r s cr T Sg) res . (3.8) 

Factorizing the scale radius r s out just translates the fact that the most important con- 
tribution to the SZ is expected to come from within a radius of r s from the cluster center 
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(therefore within a distance of 2r s along the line of sight). This is essentially true for an- 
nihilating DM, but the equation is anyway made exact for all cases by defining the average 
electron density as 



(3.9) 



This ensures Eq. fl3.8|) to be an exact formulation of the angular averaging of the optical 
depth, provided a proper definition of the dimensionless parameter l 7 n ' 7 , which carries the 
angular average: 



J 



«,7 



1 



2r s (l-/x 



res j j fji rcs 

1 ^ 



, dfi2 
2r s (1 — /i res ) j 



dfi I dl (/ 7 (r(/,/i))) r 



(3.10) 



ds (Ur(s,b))Y 



b 2 



dbb 



1 



ds / , 

y / 7 (V^T^) 



TP 



b = Dsm(9 ies ) is the impact parameter, ^ res = cos(# res ), f~ ( (r) = p y (r)/p s , s = \/r 2 — b 2 
is the reduced line-of-sight variable and R v [ r is the virial radius of the cluster. We have 
expressed the angular integral in terms of an integral over the impact parameter, which 
allows a better control of numerical aspects. It is usually convenient to integrate over b and 
s with logarithmic steps, and to define a numerical cut r cut for / 7 when b — > 0. Although 
formally set in a pure DM halo by equating the gravitational infall timescale with the 
annihilation timescale ^9|, r cut can be chosen as to ensure a good compromise between 
numerical accuracy, convergence and computing time, but should anyway obey r cut <C b res . 
We remind, nevertheless, that small core radii could instead arise from other non-trivial 



dynamical effects involving baryons, and are actually observed at the galactic scale [3C]. 

With the previous definitions, simple estimates of the optical depths associated with 
DM decay and annihilation read, respectively: 



/ 1,7 \ dec 

\ e,x / res 



/ 2,7\ann 
\ ' e,x I res 



dec ^0 
T 0,X 



T 



ann 

r °>x 777 



10 10 3 lO-^s- 1 10 17 s 

Nn 



p s /(0.05 GeV/cm 3 
m x /GeV 



10 3 kpc 



(3.11) 



{<rv) r loss 


> s /(0.05 GeV/cm 3 )" 


2 


~2r s J 2 >~i~ 


10 3 3 • lO-^cmSs" 1 10 17 s 


m x /GeV 




10 5 kpc 



where the typical optical depths are given by 



dec 

r o,x 



1.01 x 10" 



and 



ann 
T 0, X 



7.68 x 10" 



(3.12) 



It is important to note that the values that we have taken for Nq and T typify mostly 
DM particles with masses above ~ 100 MeV, since for lighter particles, the annihilation or 
decay will be mostly into light fermion pairs. For direct annihilation or decay in e + e~, we 
would have iVo = 2 and T ~ 1 (see Sec. |j), which would reduce the previous estimates of 
the optical depths by a factor of 5 x 10 3 . 
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In Fig. |], we plot the full results. We have sketched the behaviors of the exact compu- 
tation of J n "' = 2r s J nri on the one hand, and of the corresponding analytical spherical / 



line-of-sight approximations derived in Sec. 3.2 on the other hand, as functions of the res- 
olution angle (or the corresponding impact parameter, which allows to deal with physical 
distances — see the top horizontal axis). These results are translated in terms of optical 
depth by means of the DM properties used in Eq. ( p. 11 ). Note that larger values of R V i T 
would not affect the results, since the denser part of the DM distribution is within r s . 
Different values of the cluster distance D would affect the scaling of horizontal axis, the 
closer the cluster the better the spatial resolution. 



The values of the DM-induced electron optical depths derived in Eq. ( 3.11 ) have to 
be compared with what is expected for the thermal electrons. Because the latter obey an 
isothermal cored profile, as made explicit in Eq. (|2.4| ), the function J' 1,0 can be used to 
infer their optical depth (we take advantage of that we have used r c = r s ). We find: 



(r) th 

\ ' / res 



7"0,th 



n 



e.th 



0.01 cur 



where the typical optical depth is found to be 



1.0 



2r c J 
10 3 kpc 



(3.13) 



T ,th = 2.05 • 10" 



(3.14) 



with our template parameters. Such a value is very close to what is actually found in 
dedicated studies, where r ~ 10~ 3 — 10~ 2 [14]. Moreover, this value remains a constant 
function of the angular resolution up to an angular size that corresponds, roughly, to the 
core radius r c ~ 400 kpc. 

If at first glance, still disregarding the spectral aspects, we want to compare the thermal 
and DM optical depths, then we can use Fig. |l| (the left vertical axis gives 2r s J na and 
the right vertical axis provides the translation in terms of optical depth for annihilating 
DM with the generic parameters used in Eq. 3.1l| ). From this figure, we can already 



emphasize that unless (av) 3> 3 • 10 cm /s, T x 3> 10 s or m x <C 1 GeV, 7 > 1 
or p s 3> 0.05GeV/cm 3 , DM decay or annihilation cannot supply a sufficient amount of 
electrons to compete with the thermal optical depth. This even holds for rather light DM 
particles in the sub-GeV mass-scale given conventional values for the other parameters. 



Such a statement is in agreement with [10], though based on more striking analytical 
calculations, but in clear disagreement with []|, |(| |7|. The detection potential of future 
experiments with angular resolutions ~ 1" seems therefore very weak except for light DM 
with m x < 1 GeV either associated with large annihilation or decay rates, or distributed 
with very cuspy profiles (7 > 1.5). This latter feature is, however, strongly disfavored by 



most of the recent cosmological N-body simulations (e.g. [31 1). 

Before tackling the full spectral analysis, we derive hereafter some approximated ana- 
lytical expressions of the line-of-sight integral J n ^ . This set of fast formulae is aimed at 
providing some quick estimates of the optical depth for any annihilating or decaying DM 
model, and for any density profile discussed below Eq. ( |3.2[ ). 
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3.2 Analytical approximations 

The relevant part of the DM distribution is that within the angular resolution of the 
telescope. The current and foreseen angular resolutions are not better than 9 ves ~ I". 
Since nearby clusters lie at distances around 50 Mpc, the spatial resolution can hardly be 
better than b 

res ^ 

D 

@res ~ 0.1 kpc (b is defined as the impact parameter). Such values are 
much smaller than the typical scale radii r s found for cluster density profiles, so we can 
approximate the density profile of Eq. ([^) as 

P 7 (r)«p s (^) 7 * / 7 (r)«(^) 7 . (3.15) 

Beside the fact that the impact parameter is, in most cases, smaller than the scale radius 
r s , this approximation is further justified from that the dominant part of the line-of-sight 
DM contribution arises within a sphere of radius r s by the cluster center (the argument 
is weaker for decaying DM). We will therefore disregard the whole spatial extent of the 
cluster characterized by i? v i r , an d focus only on the region within r s . We further assume 
D S> r s 3> 6 rcs 3> r cut . This latter cut-off radius r cut has already been discussed below 
Eq. fl3.10| ). It has no impact except for diverging injection rates, scaling like r~ a with a > 3 
(corresponding for instance to annihilating DM distributed with a Moore profile); we will 
use r cut = 10~ 3 pc in the following. We will develop two types of approximations, one 
based on a spherical average, the other based on a more accurate line-of-sight treatment. 

3.2.1 The spherical approximation 

If the impact parameter is sufficiently small, the average value of the electron density 
determined in Eq. ( |3.9| ) can be approximated by a spherical integral of n e ^ x around 6 res , 



weighted by the conic volume ~ Tr&^es -^-vir/3 of the cluster roughly carried by the telescope 
resolution. In this case, we can derive a spherical approximation of the dimensionless 
function J" n '" f as follows 

r ra,7 ^ n-n~f ^ 12 -Rn,/3,7 ^ TO7 J 3-n 7 l^ut ' n 7 / 3 



J ~ J ^~bl s R vir r s Ts |ln(^) otherwise " 1 j 

A given product of rry corresponds to one or two DM configurations, as recalled in Tab. |2|. 
The scaling with 6 res is found to match exactly the actual behavior of l 7 n ' 7 , but we need 
to define and tune a new scale, R n ^ n , to get closer to the actual value of the optical depth 
(/3 = 3 for the NFW and Moore profiles, and (3 = 2 for cored isothermal profiles). R n R^ 



should be of the order r s for consistency reasons (see Eqs. ^8| and |3.9|) . We actually find 
that 

Rn^ = r° n J^) (3-17) 

\ 'n,7 / 

with 

r° 7 = r s ; n i7 = y and r° j7 = r 2>7 = r s , (3.18) 
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provides a reasonable approximation of the accurate numerical results, by less than one 
order of magnitude for cluster parameters not too far from those of our template model. 
This is therefore consistent with the values expected for Rn,p,>y The comparison with the 
accurate result is shown in Fig. |]. Although this approximation allows to derive orders 
of magnitude quite quickly, we recall that it is very simplistic, so one should use it with 
care. We recommend not to use it for very detailed analyses, but instead to make the full 
computation of Eq. fl3.8Q . The line-of-sight approximation that we propose in the next 
paragraph turns out to be much more accurate; still, the spherical case is rather intuitive 
and provides complementary insights. 

3.2.2 The line-of-sight approximation 



With the same approximation as above (see Eq. |3.15 ), one can also try to find analytical 



solutions to the line-of-sight integral of Eq. ( |3.10| ). Assuming that the impact parameter 
b <C r s <C<C D, which is always the case for typical angular resolutions, we can write 
J na ~ XI, where: 



'los ' 

The upper bound of the integral over s is r s and not R vu to be consistent with the validity 
domain of our current approximations. The different possible combinations of wy are 
recalled in Tab. @. 



nature / profile cored NFW Moore 

(7 = 0) ( 7 = 1) (7 = 3/2) 

decaying (n = 1) 1 3/2 

annihilating (n = 2) 2 3 

Table 2: Various combinations of the product nj and associated DM configurations for the ana- 
lytical functions an d J\^2 ( see Sees. 3.2.1 and 3.2.2). 



The function jT los can be computed analytically in the following cases: 



! (1 + jUres) . / 2 T. 

2 _ (1 + /^ros)^ r s 



res 



2 b 



JjL = (i + ih«)lr L ) in 



res 

2 



br CS J V(l + v / 2)r cu t 

Armed with these equations, the approximated optical depths are merely proportional to 
2r sl 7 1 ™J, as given in Eqs. (|3.8| ) and ( |3.9[) . The results are also reported in Fig. |l[ where 
the agreement with the exact calculation is shown to be quite good. 
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Figure 1: Line-of-sight integral for the accurate numerical computation, the spherical approxi- 
mation and the so-called line-of-sight approximation, as a function of the angular resolution over 
which it is averaged. The bottom horizontal axis is the angular resolution, the top vertical one 
is the corresponding impact parameter 6 rcs = Dsm(9 lca ). The left vertical axis is 2r s J nn — see 



Sec. 3.1 — while the right vertical axis translates those values in terms of optical depth for typical 
DM parameters, but only for annihilating DM (the optical depth axis for decaying DM would be a 
factor of ~ 1.5 smaller, for nominal parameters — see Eq. 3.11| ). 



3.3 Thermal to DM optical depth ratio i] x 

To size the relative amplitude of the SZ distortion induced by DM annihilation or decay on 
top of the thermal component, it is useful to derive the ratio r\ x = (T ejX ) ros /(r eit h)res ^ 0.1, 
at any frequency where AJ t h ^ 0. Since the frequency where AI crosses is different 
for each electron component, it is further important to compare the absolute amplitude 
of A/ x , at the frequency where A/ t h = 0, with the current experimental sensitivities, i.e. 
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AI/I ~ 1CT 6 . We will estimate the ratio r\ x below, but postpone our discussion of the 
latter point in Sec. ||. For the sake of clarity, we will base the following calculation on the 
assumption that the DM density profile is an NFW (7 = 1). Likewise, we will use the 
analytical line-of-sight approximation derived in the previous section, which turns out to 
be the most precise. 

Using Eqs. ( |3.11| ) and ( |3. 13 ), the ratio of the DM to thermal optical depths is propor- 



tional to J\qs' / 3\os-> f° r decaying (annihilating) DM. More explicitly, we get, for decaying 
DM, 



4 ec 



I _l,l\dec 
\'e,X/res 

(''"e,th)res 



.-dec 

r o,x 



TQ.th 



n 



e,th 



0.01 cm- 



iVo T_ r x 

10 10 3 lO-^s- 1 10 17 s 



71o 



^O.dec 



-1 -1 



n 



e,th 



In 



0.01cm~ 3 

2(r s /400kpc) 
6 res /lkpc 



Tlo 



^0 t F * 

10 10 3 lO-^s" 1 10 17 s 



p s /(0.05 GeV/cm 3 ) 
m x /GeV 

Ps/(0.05 GeV/cm 3 ; 
m x /GeV 



J} 



los 



7° 



(3.21) 



with 



r/°' dec = 3.29 x 10" 4 



For annihilating DM, the result is proportional to J\ os /J\ os , and we have instead 



„ann 



/ l,l\ann 
\'e.X/res 

(7"e,th)res 



0,ann 
' v 



e,th 



0.01 cm~ 3 
r s /400kpc' 



(av) 



Tlo 



10 10 3 3 • lO^cmSs- 1 10 17 s 



p s /(0.05 GeV/cm 3 
m x /GeV 



frres/lkpc 



where 



0,ann 
''X 



4.71 x 10' 



(3.23) 



(3.24) 



The relevant parameters here are ?7 x dec and ?y x ' ann , which both are <0.1 for most of 
DM models and for current and future experimental performances. Because the optical 
depth sizes the amplitude of the spectral distortion (see Eq. |2.1| ), this means that detecting 
a SZ signal from DM annihilation or decay on top of the thermal contribution demands 
very strong spectral distortions of the former with respect to the latter. We study the 
spectral features of the SZ signatures in the next section. 



4. Spectral distortion analysis 



To complete our study, we have to tackle a full spectral analysis. Nevertheless, what is 
quantitatively interesting is the amplitude of the spectral change due to DM at the SZ 



- 11 - 



transition frequency of the thermal gas, i.e. when its global effect is null. In this section, 
therefore, we estimate the amplitude of AI X at the frequency for which A/ t h = 0, to size 
the potential impact of DM at the transition frequency of the thermal background. We 



adopt the formalism developed in [32], which is very well suited for our calculation. We 
refer the reader to that article for more details on the quantities and results that are used 
below. 

Let us first define the reduced blackbody intensity i(x) = I^E^/iq = x 3 /(e x — 1), 
where x = E^/^kTo), Tq is the CMB temperature, and iq = 2(/cTo) 3 / (he) 2 . With these 
notations, Eq. (|2.1|) can be rewritten in terms of the reduced intensities 



5i(x) = r(j(x) -i(x)) , (4.1) 
where we further define the scattered photon intensity as: 

Z 1 ^- m ax 

j(x) = / dtP!(t)i(x/t) . (4.2) 







The photon frequency redistribution function P\(t), valid in the single scattering ap- 
proximation, depends on the normalized injected electron spectrum J-(E) = J-(p) (p = 
p/(m e c) = 7e/3 e = \/le ~ 1 being the reduced electron momentum) as follows: 

Pi(t) = J dpFffiPtfrp), (4.3) 

where P\(t,p) has the following analytical form 

Pl{t,P) = ~~32pW i 1 + ( 10 + 8 P 2 + ^)t + t 2 }+ (4-4) 
3(1 + t) [ 3 + 3p 2 +p 4 3 + 2^/ \Ht)\\\ 

The maximal frequency shift, above which Pi(t,p) = 0, fulfills the condition |ln(t max )| = 
2 asinh(p). 

Armed with these equations, we can derive the intensity shift 5i for any electron com- 
ponent, given its spectrum. For the thermal electrons, we will use a Maxwell-Boltzmann 
spectrum 

•FthGO = J^jP 2 ex P (-/Wl+P 2 ) , (4-5) 

where (3th = m e c 2 / '(fcTth) is the inverse reduced electron temperature, and where the nor- 
malization is ensured with the modified Bessel function of the second kind Ki . The resulting 
redistribution function (see Eq. |4.3j ) is shown in Fig. 0, for different temperatures. 

For the DM electron yield, we will assume that the injected spectrum is oc 5{E — 
nm x /2), with n = 1 or 2 for decay or annihilation, respectively. This is equivalent to 
assuming a direct annihilation or decay in e + e~, which is particularly relevant for light 
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Figure 2: Frequency redistribution function for different thermal electron gas temperatures, fclth 
5, 10, 15 keV. 



DM particles but a very optimistic spectral hypothesis in the general case. We can thus 



determine the normalized equilibrium spectrum with Eq. (3.5): 



HE) 



K 

W) 



(for E < nm x /2) , 



(4.6) 



where K = K(E m i n , E ma ^) is a normalization constant such that J dET{E) = 1 in the 
energy range of interest. For the energy loss rate b(E), we can combine the two most 
relevant regimes, i.e. the Coulomb losses below Eq = 1 GeV, of timescale rc ui. an d the 
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Figure 3: 5ith(x)/r for different thermal electron gas temperatures, kT t h = 5, 10, 15 keV (the right 
vertical axis gives |(Jith(^)| for our template cluster, with an optical depth of r = 2.05 x 10~ 2 ). 



inverse Compton losses on CMB above [33], of timescale tic- Actually, we have: 



b(E) 



10- 16 GeV/s 



bi C {E) + bconl(E) 



E J ft 1cm J 



(4.7) 



1 + ln 



E/E 



n e ,th/lcm 3 



where the Coulomb loss rate is taken from [34|. We can simplify the previous expression 
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by means of the timescales 

TCoul. & Tic r Coul. [ \^0/ j 

T Coui. _ n 1 a w /3e tic _ o w i _ r coui. /< \ 

— = 0.16 x ; 5- ; — = 3.77 and a T = . (4.8) 

10 16 s n ejth /lcm- 3 ' 10 16 s nc 

With this form of energy loss rate, we can reexpress Eq. (|4.6j) as 



F{E) = K ( E »^>-x> . (4.9) 
a T #7 + 1 

The correspondence between energy and reduced momentum is straightforward: p = 
1/7I — 1, and dE = pm^c dp/E. The energy lower bound E m i n must be consistent with 
the cluster age, DM-induced electrons having essentially been injected since the cluster 
has formed. Therefore, these electrons had only a time of the order of the cluster age 
At ~ t c \ to lose their energy. We can compute the minimal energy E m [ n > m e c 2 corre- 
sponding to any injected energy E by demanding dE' < At, which is equivalent 
to K^ 1 (E m i n , E) < At. For typical cluster formation redshifts around 5, E m - m is merely 
found to be m e c 2 . 

The photon redistribution function obtained from the DM-induced electrons, once the 
normalized spectrum defined in Eq. ( [4.6D has been injected into Eq. (|4.3|), is shown in 
Fig. [|, for different WIMP masses, assuming direct annihilation or decay in e + e~ . 

Notice that when coming to the numerical calculation of the optical depth, we will 
take advantage of the estimates performed in the previous section by remarking that the 



function T defined in Eq. (3/7) obeys rigorously Ti oss .F = K~ l {E m \ n , m x ), in the current 
DM configuration. If we assume that ti oss ~ t cou i ~ T{ c , then we have T ~ 1. Hence, with 
such an assumption, we can rescale the optical depths given in Eq. (|3.12| ) accordingly. 

Our aim is to size <5i x (xth) at a frequency xth such that <5ith(^th) = 0. x t h is readily 
found by computing Eq. (f4.1| ), given the thermal electron temperature. For completeness, 
we take three different temperatures, kT t h = 5, 10, 15 keV, though it is well known that 
the transition frequency is much less dependent on the temperature than the amplitude, 
as shown in Fig. |3[ As illustrated in that plot for the three cases, we find x t h — 4. We can 
now calculate 

5i x ( x th) = T x (j x (x th ) - i(x th )) , (4.10) 

with 



j x (xth)= dtPy(t,p x )i{x/t) , (4.11) 







and for different DM particle masses. 

We first report the behavior obtained for 5i x (x)/r x with different WIMP masses in 
Fig. [|, assuming direct annihilation in e + e~ (also valid for decaying DM, but corresponding 
then to twice heavier WIMP masses). The factor of 1/t x allows a prediction independent 



- 15 - 



1 — 



10" 1 



io 2 — m v = 100 MeV 



10 ; 



^10* 



10' 



10 1 



10 7 



10' 



10 



iiiiiiiii iiiiiiiii 1 1 iiiiii iiiiiiiii 1 1 iiiiii 1 1 iiiiii iiiiiiiii iiiiiiiii iiiiiiiii iiiiiiiii lining iiiiiiiii iiiiiiiii iiiiiiiii iiiiiiiii inn 



Redistribution function 



m 7 = 1 MeV 



m y = 10 GeV 




I I I ml I I I M I 1 1 1 1 1 1 1 1 1 1 



10" 10' 7 10" 6 10" 5 10" 4 10" 3 10* 10' 1 1 10 1 2 1 3 1 4 1 5 1 6 1 7 



Figure 4: Redistribution functions for different WIMP masses, assuming \X e+e ■ Decaying 
DM functions are identical for \ e + e~ but with twice the indicated masses. 



of the electron density, i.e. independent of the cluster halo parameters. The right vertical 
axis gives the conversion in terms of 5i x (x), assuming t x = 1.5 x 10~ 9 (see below), which 
reaches a maximum of 5i ~ 10~ 10 around the transition frequency of the thermal component 
xth ~ 4. This maximal value corresponds to an annihilating DM of 1 GeV distributed with 
an NFW profile and observed from our template cluster within an angular resolution of 
10" (see Eq. |T| and Fig. |l]). It is anyway easily rescaled for other configurations. 

In Fig. ^, we also trace <^ x (x t h)/(T x m™io(x t h)), i-e. at the transition frequency of the 
thermal component, as a function of the DM mass m x ; n =1 (2) for decaying (annihilating) 
DM. We implement the expected mass dependence of the optical depth thanks to the factor 
of 1/m™. To compare with experimental sensitivities, we have converted our result in terms 
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of the relative intensity shift 5i, but only in the case of DM annihilation (cf. the right vertical 
axis). To do so, we have used the angular average of t x as derived in Eq. ( |3.11 ), which rests 
on the template DM halo defined in Tab. [|. Likewise, to be consistent with the assumption 
of annihilation or decay in e + e~, however, we have used iVo = 2, and, as discussed above, 
T = 1 (i.e. for to iX , we have taken the values of Eq. 3.12 divided by 5 • 10 3 , which gives 
To iX = 1.5 x 10~ 9 ). To switch to the Si associated with DM decay, one may apply a factor 
of 1.01/7.68 ~ 0.13 to the right vertical scale (see Eq. 3.11| ). Notice that the scaling with 
m x is trivial, oc l/m x 2 for DM annihilation, and oc l/m x for DM decay. 

From the results presented in Figs. [5] and ||, it seems quite difficult to detect any SZ 
signal from DM annihilation or decay, since the typical values obtained for the relative 
intensity shift is 5i x ~ 10~ 10 , still far from current experimental sensitivities (Planck was 
optimized for 5i ~ 10 -6 [11], i.e. the /uK level in terms of temperature). However, we 
remind that this particular value of the intensity shift is connected to some assumptions. 
For the line-of-sight integral, which depends on the cluster halo properties, we have used 
1r s J = 10 3 and 10 5 kpc for DM decay and annihilation, respectively, which corresponds 
to observing our NFW template cluster within an angular resolution of ~ 10" (see Fig. [I]). 
Moreover, there are additional assumptions on m x ~ 1 GeV and on the DM properties, 
though quite generic (see Eq. 3.11 ). Anyway, we emphasize that these assumptions are 
rather optimistic as regards the DM modeling, and should therefore be considered as con- 
servative. For comparison, the intensity shift is shown to be much larger for thermal 
electrons (see Fig. || with the vertical right axis), 5i ~ 10 -4 — 10 -3 , over a broad band of 
frequencies, except at the transition frequency. Hence, the DM effect is likely too small in 
the major part of the "natural" DM parameter space, but might be still sizable for very 
light DM candidates and very cuspy profiles. 

To figure out the influence of both the DM halo cuspiness and the experimental angular 
resolution, one can read off the left vertical axis of Fig. [l], which provides values of the line- 
of-sight integral for different cluster halo profiles as functions of the angular resolution. To 
boost our predictions by > 5 orders of magnitude we would have to consider very cuspy 
profiles (7 > 1.5) and in the meantime very deep spatial resolutions (6 rcs ~ 0.1 kpc, or 
equivalently 9 res ~ 1"), at least for the template set of halo parameters of Tab. |]. We will 
show in Sec. ||, however, that even this over-optimistic case is actually limited by spatial 
diffusion effects. 

If we stick to an NFW profile, an alternative to enhance our predictions would be 
to consider very light annihilating WIMPs, in the sub-GeV mass-scale. However, in that 
case, the value of the annihilation cross considered in Eq. ( |3.11 ) is bounded by astrophysical 
observations in the Milky Way. Indeed, it was found in |35|, [36| , that the annihilation cross 
section for light DM must obey (av) < 10~ 31 (m x /MeV) 2 cm 3 s _1 to not overshoot the 511 
keV bulge emission measured by the SPI spectrometer onboard the INTEGRAL satellite 
1 37 ] . Since Fig. || was done with a cross section of (av) = 3 • 10 _26 cm 3 s~ 1 , such a limit 
translates into a plateau below m x ~ 500 MeV, translating into a saturation of 5i x < 10 -8 , 
much below the current and future sensitivities. 

Invoking very cuspy profile with 7 > 1 appears be the only possibility left for DM to 
generate an observable SZ distortion of the CMB spectrum, but such an extreme configu- 
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Figure 5: \5i x (x)\/r x = \j x (x) — i(x)\ as a function of the reduced frequency x for different WIMP 
masses, assuming direct annihilation in e + e~ (identical for decaying DM, but corresponding to 
twice the masses). The right vertical axis provides the conversion to |<5i x (x)| for an optical depth 
of t = 1.5 x 10 -9 , corresponding to a configuration with 1 GeV annihilating DM with an NFW 
profile inside our template cluster observed with an angular resolution of 10" — see Eq. ( 3.11 ) and 
Fig. 0. 



ration (i) is hardly motivated from theoretical and observational constraints and (ii) will 
even be shown, in the next section, to have limited effect due to spatial diffusion. Last 
but not least, one should not forget about the additional foreground coming from other 
relativistic electrons also injected in clusters from standard astrophysical sources [32]. In 
any case, any attempt of DM interpretation of any SZ imprints would seem daring, at least 
in our view. 
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Figure 6: 5i x (x t h) / {m™T x ) as a function of the DM particle mass m x , for x t h = 4 (left vertical axis). 
The right vertical axis provides the conversion to |5z x (xth)| for an optical depth of r = 1.5 x 10~ 9 , 
corresponding to a configuration with annihilating DM with an NFW profile inside our template 
cluster observed with an angular resolution of 10" — see Eq. (fTO) and Fig. EL 



5. Impact of spatial diffusion: analytical insights 

So far, we have assumed that the equilibrium electron density was set only from energy 
losses and neglected the spatial transport, which was presented as a reasonable hypothesis 
in M for angular resolutions < V . Nevertheless, as discussed above, angular resolutions 
down to 1", as expected in future experiments, correspond to quite small spatial scales for 
typical nearby clusters, so that the relevance of such an assumption can be questioned. 

If we include spatial diffusion, then Eq. ( |3.3| ) must be rewritten in terms of a current 



conservation equation 



-V 



dn \ 
dEJ 



Q(x,E) . 



(5.1) 



K r is the diffusion coefficient that characterizes the stochastic transport of electrons caused 
by their diffusion on the magnetic inhomogeneities of the intracluster medium. We will 
assume this coefficient to depend on the energy only K r = K r {E). This equation can be 
solved by looking for the associated Green function, which is readily found in infinite 3D 
space: 



Q(x,E «- x s ,E s 



1 



■ cxp 



X 



b(E)TT2\3 

where A is a characteristic propagation scale defined as 



A 2 



(5.2) 



X 2 (E,E S 



Je b(E') 



(5.3) 



The electron density at any position x and energy E in the cluster is then given by 



dn 
HE 



(x,E) 



n m x /2 



dE s I d 3 x s G(x, E «- x s , E s ) Q(x s ,E s 



(5.4) 



We can remark that we recover the diffusion-less case in the limit of vanishingly small 
propagation scale A — > 0, since in this case we have G(E, x x s ) 5 3 (x — x s )/b(E). For 
a monochromatic injection of electrons from DM annihilation or decay, the diffusion- less 
limit is therefore recovered when E — > nm x /2. More generally however, since the Green 
function exhibits a Gaussian behavior of typical scale A, the electron density is expected 
to be smeared out over this scale, which we actually demonstrate below. 

The Green function appearing in Eq. ( |5.2| ) was derived assuming diffusion in an infinite 
3D space, which is obviously not the case since clusters have finite sizes. This still holds 
provided the propagation scale A is much lower than the typical size of the cluster, and 
while electron injection at the border of the object is irrelevant compared to injection in 
the more central regions. If set by magnetic inhomogeneities, the diffusion coefficient can 
be expressed as a function of the regular component of the magnetic field. Although it is 
absolutely not clear whether or not a Kolmogorov spectrum is relevant to describe magnetic 
turbulence in clusters, since it is not the case at the Galactic scale, it is still conventional 



to adopt such a behavior for the diffusion coefficient [39]: 
K r (E) = K 



\ Ib ' 


2/3 


" B ' 


-1/3 


E 


1/3 


20kpc 








GeV 





(5.5) 



where B is the magnetic field and Ib its coherence length. With the above values, which 
we will use in the following, the normalization of the diffusion coefficient can set to Kq ~ 
2.3 x 10 29 cm 2 s _1 [39]. Since the typical energy loss timescale is of order ti oss ~ 10 16 s, 

' 30kpc <C r s , which is actually quite large though still smaller 



we have A ~ t\ 
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than the typical scale radii of clusters — note that this diffusion scale will slightly decrease 
with energy above 1 GeV, when the inverse Compton dominates over the Coulomb losses. 
Therefore, because the injection rate of electrons is expected to dominate within a volume 
set by r s <C -R v ir and since A < r s , we can safely disregard the radial boundary condition 
and use the infinite 3D solution written above. 

The source term Q is defined in Eq. (|3.1|). To simplify the discussion, we will further 
suppose that DM annihilates or decays into electron-positron pairs, so that Nq = 2 and 
F(E S ) = 5{E S — nm x /2) — n = 2 in the case of annihilation, 1 in the case of decay. In this 
fiducial instance, the propagation scale only depends on the energy E, A = X(E, n m x /2) - 
see Eq. ( |5.3D . Therefore, the electron density at any cluster radius r and energy E < nm x /2 
reduces to 



dn 
dE 



r,E) 



4 a n p r s 



b{E)^X* J_ x 



dp! I dr' r 



cxp 



r' 2 + r 2 — 2 r r' p! 



A 2 



/>'), (5-6) 



where we have taken advantage of the spherical symmetry of the source term, the radial 



dependence f-y(r') of which has been defined in Eq. (3.2). The integral over the cosine p! 
is straightforward, and we are left with: 



dn , „. 

HE {r < E) 



2 a n p r s 



b{E)^Xr 

2a nP n s rT 
b{E)^Xr 



k=l 

£(- 

k=0 
k=l 

E(- 

k=0 



dr' r' f~(r') exp 



(r 



-l) k r') 2 



A 2 



(5.7) 



1 



r+A 



dr'\r'\ l - n ^ exp 



l) k r') 2 



A 2 



In the last line, we have employed the approximation given in Eq. (3.15), where the DM 
profile is taken as a simple power law; indeed, we are mostly interested, again, in the SZ 
contribution from the densest central parts of the cluster. Moreover, we have restricted 
the integral range to [r — A;r + A], accounting for the fact that electrons coming from 



distances above A are Gaussianly depleted, consistently with Eq. (5J2); we remind that 
A ~ 10 kpc. In the very central region of the cluster and within this range for r', the 



arguments in the exponentials 



■1) £ A > so that 

we can expand the integral to 



get more insights on the physics at stake. To the first order, the result is analytical, and 
we find after integration: 



dE { ^ E) 



(r + A) 



3— n-y 



1 — sign(r — A) 



3— n-y 



(r + A) 



3— n-y 



b(E)^X 3 (3-n 7 ) 

where we let the reader derive the logarithmic expression that arises when 3 — 717 = 0. This 
result is very interesting, since it states that at the center of the cluster, i.e. r = 0, the 
electron density is no longer diverging like r - ™ 7 when diffusion occurs, but instead behaves 
like: 

WanpU^-nj)- 1 -■ 



b(E) 



TT 



X 



n-y 



(5.8) 



As expected, the electron density scales like A n7 , where A is the diffusion scale, but we have 
grounded this intuition on a rather rigorous demonstration. Of course, this approximation 
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is not valid for vanishingly small values of A. This means that at any energy but the 
maximal energy injected from DM annihilation or decay — when A is sizable — the actual 
electron density is not governed by the (squared) DM density itself, but is diluted away 
over a region delineated by A. 

This has important consequences for the SZ calculation. Indeed, we have computed the 
optical depth in the diffusion-less approximation, showing that it was significantly increas- 
ing with angular resolution for cuspy profiles. Nevertheless, we have just demonstrated 
that the electron density should in fact saturate within a scale fixed by ~ A, leading to the 
saturation of the optical depth itself, independently of the angular resolution. For A ~ 10 
kpc, this would correspond to an angular size of ~ 40" for our template cluster, as depicted 
in Fig. H|, which means that instruments with better resolution than this value would not 
improve the detection potential. 

The full SZ computation is slightly different when spatial diffusion is considered, since 
we can no longer factorize the energy and the spatial parts as we did in Eqs. (^5,3.5,3^6|), 
but must use Eq. fl5.4|) to compute the the optical depth. Nevertheless, up to a good 
approximation valid in the monochromatic injection case, we can associate an energy- 
dependent core radius of size \(E,nm x /2) to each halo model. Given the values that we 
used for the diffusion coefficient and for the energy loss rate, this would lead to a plateau 
in the optical depth curves for scales below ~ 10 kpc, as also found in the final published 
version of [ 10 1 . Since our spectral analysis was performed for an NFW profile and an 
angular resolution of 10", the full spatial diffusion calculation would not only decrease 
the predicted amplitude by a factor of ~ 2, as seen from Fig. |], but also make this result 
constant for any better angular resolution. Should have we done this exercise with a Moore 
profile, we would have overshot the actual amplitude by one order of magnitude. 



6. Conclusion 



In this paper, we have revisited the prospect for radio experiments, like Planck or ALMA, 
for observing SZ distortions of the CMB spectrum due to the presence of relativistic elec- 
trons induced by DM annihilation or decay. Whenever possible, we have provided analytical 
grounds to the discussion. 

In Sec. H|, we have first focused on the most important physical weight on the distortion 
amplitude, the optical depth associated with the electron density (the thermal background 
as well as electrons induced by DM annihilation or decay). We have considered typical 
values for the DM particle properties and a generic template cluster located at a distance 
of 50 Mpc (Coma lies at ~ 90 Mpc), with features very close to those observed in N- 
body simulations and a thermal electron population with a density of n eit h = 0.01cm -3 , 
consistent with most of X-ray observations. 

In the frame of the diffusion- less limit proposed in p], we have stressed the importance 
of the angular resolution inside which the SZ signal has to be averaged, showing that even 
for very optimistic values of ~ 1", the thermal electron optical depth almost always strongly 
dominate over that of the DM-induced electrons. This actually means that the thermal SZ 
completely overcomes the DM SZ at all frequencies but the thermal transition frequency. 
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In addressing the angular resolution issue, we have derived simple and useful analytical 
approximations of the line-of-sight integral (Eqs. 3.16 and 3.20| , the latter being more 
accurate) that can be used to compute the optical depth for any decaying or annihilating 
DM model. Actually, they are also useful for quick predictions of 7-ray fluxes from DM 
annihilation or decay in any extragalactic objects, subhalos, dwarf spheroidals, galaxies or 
clusters. 

Then, we have performed the full spectral analysis of the SZ distortion due to DM- 
induced electrons in Sec. [|, showing that it can hardly exceed the level in terms of 
temperature fluctuations, even for e.g. very light DM particles, which are further con- 
strained by 511 keV observations of the Galactic center. Unconventional very cuspy cluster 
halo profiles and very large annihilation cross sections or decay rates that would escape 
other astrophysical contraints for some reasons might still lead to signal close to future 
sensitivities. However, such configurations are not supported by current studies in each of 
the specific inputs and are therefore quite unlikely. 

Finally, we have discussed the impact of considering the spatial transport of electrons 
in clusters in Sec. ||, providing again analytical results. We have demonstrated that the 
diffusion-less limit is actually not valid in the cluster center for SZ computations with small 
angular resolutions because the electron density is actually smoothed over a scale set by 
the transport scale. This scale is small only for energies close to their injection values, 
and are therefore sizable over most the whole energy range but for E — > m x in the case of 
annihilating DM. Given the propagation parameters used in our analysis, this scale is of 
order ~ 10 kpc, below which the optical depth saturates. This means that improving the 
angular resolution below that scale will not improve the detection potential, unfortunately. 
Our results are in good agreement with those obtained in [1C], to which they provide deeper 
analytical grounds. Reversely, we cannot support the statements made in [||, |(| f?J about 
the promising potential of SZ observations for indirect detection of DM. 
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